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We introduce a class of higher-order derivative models in (2, 1) space-time dimensions. The models 
are described by a vector field, and contain a Proca-like mass term which prevents gauge invariance. 
We use the gauge embedding procedure to generate another class of higher-order derivative models, 
gauge-invariant and dual to the former class. We show that the results are valid in arbitrary (d, 1) 
space-time dimensions when one discards the Chern-Simons and Chern-Simons-like terms. We also 
investigate duality at the quantum level, and we show that it is preserved in the quantum scenario. 
Other results include investigations concerning the gauge embedding approach when the vector field 
couples with fermionic matter, and when one adds nonlinearity. 
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O ' I. INTRODUCTION 

(N 

It is hardly necessary to remind the remarkable and powerful properties of the duality mapping as an analytical 
^ ■ tool in Field Theory as well as in String Theory 0]. On the other hand, the interest in the study of theories involving 
^ higher-order derivative is by now well-appreciated and remains intense. Within the context of the Maxwell theory, 
generalizations involving higher-order derivatives can be found in More recently, in Refs. one finds 

generalizations that involve both the Maxwell and Chern-Simons (CS) terms 0,0, 0. In the present work we study 
. duality symmetry in extended theories, which contain higher-order derivatives involving both the Maxwell and CS 
K>" ' terms. 

The interest in the subject has been recently fueled by several motivations, in particular by issues related to String 
Theory. As one knows, string theories engender the feature of containing interactions that may involve an infinite 
number of space-time derivatives, and that may lead to standard Field Theory in the low energy limit. Thus, higher- 
order derivative contributions would certainly appear when one consider next-to-leading orders in the energy. Based 
^sj I upon such interestin g po ssibilities, there is renewed interest in investigating higher-order derivative models, inclu ding 
, the case of tachvons |lOl HTl IT^ llsL fl^ ITsl ITgL IT^ . Other lines of investigations can be found for instance in T^, 1^ 
and also in gravity, where higher orders in the scalar curvature R are considered in the so-called nonlinear gravity 
theories - see, for instance, the recent works |23, 0| and references therein. 

• The presence of higher-order derivatives to control the behavior of systems is not peculiar to String and Field 
, theory. It may also appear in condensed matter, and can be important to describe higher-order phase transitions 

• |23, as for instance the case of the fourth-order transition in superconducting Bao.6Ko.4Bi03, which is described in 
[23|. Furthermore, sometimes one has to include higher-order derivative of the order parameter to correctly describe 
pattern formation in chemical reactions, and in other branches of nonlinear science 24, 2^. 

We investigate the subject dealing with issues that appear very naturally in Field Theory. Specifically, we examine 
the duality mapping of higher-derivative extensions of self-dual (SD) and Maxwell-Chern-Simons (MCS) theories, 
including the presence of fermionic matter. We work in the (2, 1) dimensional space-time, with e*^^^ — £912 = 1; our 
metric tensor has signature (-I-, — , — ) and we use natural units. We start in the next Sec. ^ where we examine the 
dual mapping of the MCS-Proca model with the MCS-Podolsky theory. Our methodology makes use of the gauge 
embedding procedure, an approach which has been shown to work very efficiently to unveil the dual partner of a 
specific model [2^ |23, |2^ 113 113 ■ The investigation of the MCS-Proca model is new, and we use it to set stage to 
generalize the model to higher-order derivatives. We develop this generalization in Sec. IIIII where we introduce the 
main model, which is defined in terms of n = 1,2, and in Sec. IIVI where we consider the case n — + 00. We show 
in Sec. IIIII that if one discards the CS and CS-like terms, the results are then valid in arbitrary (d, 1) space-time 
dimensions. We also show in Sec. IIVI that the duality is preserved at the quantum level. In Sec. we examine the 
presence of matter, coupling fermions to the system. In Sec. lVII we change the model introduced in Sec. 11111 to include 
the presence of nonlinearity. We end our work in Sec. IVIII where we present our comments and conclusions. 
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II. DUALITY TRANSFORMATION IN THE MCS-PROCA MODEL 



It is well know that both the SD and the MCS models are dual representations of the same dynamics: they carry 
one massive degree of freedom of definite helicity determined by the relative sign of CS term. However, the SD 
representation hides a gauge symmetry, which is explicit in the MCS model. This is easily seen when one establishes 
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the correspondence ^ ^ e^^pd^AP, which maps the SD field into the dual of the basic field of the 
MCS model. We notice that in the above dual mapping, one relates a gauge non-invariant model with an equivalent, 
gauge invariant theory. In this process, however, the non-gauge field is identified with a special form containing the 
derivative of the gauge field of the dual model. The identification of this mechanism is crucial for the generalizations 
that we will implement below. 

In this section we elaborate on a prototype of the study we intend develop in this paper. We discuss the theoretical 
motivations for studying higher derivatives dualities, present a physical scenario for the applications of the ideas 
elaborated and the quantum implications. 

In the past, there have appeared several different ways of extending electrodynamicSjtrying to smooth infrared or 
ultraviolet singularities that appear at large or short distances. One is the Born-Infeld type of generalization, and 
others include the generalizations introduced by Proca Q and by Podolsky The Born-Infeld approach involves 
nonlinear extension, which will be considered in Sec. IVII The Proca model involves the addition of a mass term for 
the vector field, which was introduced to smooth infrared singularities. The Podolsky model involves higher-order 
derivatives and was introduced to smooth ultraviolet singularities. Thus, the Proca and Podolsky model deal with 
dual aspects of the electromagnetic interaction, and so they may be connected by some duality procedure. 

The dual aspects of electrodynamics that appears in the Proca and Podolsky models constitute the central subject 
of the present Section. This duality is also of importance in the extension of the bosonization program from D=2 to 
higher dimensions. The mechanisms of the 3D bosonization in particular are quite dependent on the duality results 
involving the presence of the Chcrn-Simons term. Besides providing us with the proper scenario for the applications 
of our ideas, bosonization will also be crucial for the interpretation of the new parameters in the Podolsky extension. 

A. The duality Procedure 

In order to exemplify the general procedure of duality, we consider the more general model 

L = —A^^A,, - -F^^F^, ~ -me^,xA^d-A>^ (1) 

where m is the mass parameter. We have introduced the real and dimensionless parameter a in order to obtain the 
SD model for a = 0, or the MCS-Proca model for a = 1. 

The equation of motion involves second-order derivatives, and the model is supposedly dual to a g eneralized MCS- 
Podolsky theory. To verify this assumption we follow the gauge embedding procedure [2a. I27L |28L l29l l30l | to construct 
its dual equivalent model. Firstly, we compute the Euler vector associated with the MCS-Proca theory. We get 

= m^A^ - m Sf^^xd^A^ + a d^F,^ (2) 

The first iteration leads to 

£1 = £0 - K^'B^ (3) 

where Co is identified with MCS-Proca Lagrangian given by Eq. Q). Also, i?^ is an auxiliary field, which varies 
according to 

6B^ = SA^ = d^A (4) 

This choice makes the non-invariant term in Cq to cancel with the term K^SB^. Therefore 

2 

6Ci = -B,,SK^ = -^-6{B^) (5) 
Thus, we can write the gauge invariant second-iterated Lagrangian 

£2 - £0 - Kf^B^ + (6) 
This ends the iteration process, and we can eliminate the auxiliary field to obtain the dual model 
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or better 




models. We recall that the Podolsky model was introduced in order to smooth ultraviolet singularities. In this sense, 
we notice that our approach for duality is working standardly, since we are linking infrared and ultraviolet problems, 
in MCS-Proca and MCS-Podolsky models. 

In the above investigation, if a = the MCS-Proca model in becomes the standard SD model; in this case, in 
the generalized MCS-Podolsky model ^ one kills the Podolsky terms, and we are led back to the MCS model. On 
the other hand if a = 1 we get to an extended CS-Podolsky model. 

It seems important, at this juncture, to establish the scenario at which the parameter a becomes a physical quantity. 
This is done next where we show that action (jSJ for the generalized MCS-Podolsky model is the low energy effective 
action for a self-interaction fermionic model. To set the problem in a proper perspective, we review shortly the 
bosonization procedure. 



It is well known by now that the bosonization in D=3 mapps a massive scalar particle coupled to a Chern-Simons 
gauge field into a massive Dirac fermion for a special value of the Chern-Simons coupling. This is a relevant issue in 
the context of transmutation of spin and statistics with interesting applications to problems both in quantum field 
theory and condensed matter physics. This boson-fermion transmutation is a property which holds only at very long 
distances, namely, at scales long compared with the Compton wavelength of the particle. As so these results hold to 
lowest order in an expansion in powers of the inverse mass of the particle. 

The equivalence of the three dimensional effective electromagnetic action of the CP^ model with a charged massive 
fermion to lowest order in inverse (fermion) mass has been proposed by Deser and Rcdlich ,31j. Using the results of 
[3l|. bosonization was extended from two to three dimensions 32j. However, contrary to the two-dimensional case 
where the fermionic determinant can be exactly computed, bosonization in higher dimensions is not exact and, in the 
general case, it has a non-local structure. However, for the large mass limit in the one-loop of perturbative evaluation, 
a local expression materializes. Indeed, it has been established [s^, to the leading order in the inverse fermionic 
mass, an identity between the partition functions for the three-dimensional Thirring model and the topologically 
massive U(l) gauge theory, whose dynamics is controlled by a Maxwell-Chern-Simons action. Here we show that the 
contribution next-to-leading order leads to the MCS-Podolsky model via duality transformation. 

Below we compute the low energy sector of a theory of massive self-interacting fermions, the massive Thirring 
Model in 2 -I- 1 dimensions, that can be bosonized into a gauge theory, the Maxwell-Chern-Simons gauge theory and 
its possible higher derivative extensions. We start from the fermionic partition function for the three-dimensional 
massive Thirring model, following closely the spirit of j32| , 



with the coupling constant having dimensions of inverse mass. Here tp^ are N two-component Dirac spinors and 
is a global U(l) current defined as, 



Notice that we have reverted to Euclidean metric in this subsection. As usual, we eliminate the quartic interaction 
by performing a (functional) Legendre transformation through the identity 



B. Physical Interpretation 




(9) 



J 



(10) 




(11) 



after a scaling gA^/^/N. The partition function then becomes 




(12) 



Formally, the fermionic path- integral gives the Dirac operator determinant. 




(13) 
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The determinant of the Dirac operator is an unbounded operator and requires regularization. Bosonization will 
depend on the actual computation of this deteminant, namely, whether it can be computed exactly in a closed form 
or an approximate recipe must be enforced. In general it leads to non-local structures but, under some approximation 
scheme (like the inverse mass), a local result emerges. 

This determinant can be computed exactly for D — 2, both for abelian and non-abelian symmetries. For the D=3 
case this determinant has been computed in 31] as an expansion in inverse powers of the fermion mass giving, in the 
leading order, the Chern-Simons parity violating term, as well as parity conserving Maxwell term, which is central to 
our discussion here, 

lndet(^^ + M+-^4^ =±^ Z'^''"-^'''^"^'^" 2^M j d^'xF'^'' ^ 0{d^ /M^) (14) 
We bring these results into the partition function to get 

ZTh= ypA^e-^^^I^-l (15) 

where Seff[A^] is given, up to order 1/to, by 

Sef}[A,] ^IJd'x (^A,A^^ T - ^F'^'^F,,^ (16) 

which, after the scaling — s- m^^, with m = ^iv/g^. With these identifications we find a — In conclusion, we 
have established the following identification: 

Zrh ~ ZMCS-Proca (17) 

which is valid to next-to-leading order in 1/M. 

In the preceeding section we have established the dynamical equivalence between the model defined by SMCS-Proca 
the MCS-Podolsky theory. This proves the equivalence, to this order in 1/M expansion of the partition functions for 
Thirring model and the MCS-Podolsky theory: 

Zrh ~ Zmcs- Podolsky (18) 

It expresses the equivalence between a theory of 3-dimensional interacting fermionic theory and Maxwell-Chern-Simons 
vectorial bosons in the long wavelength aproximation. It is interesting to observe that the Thirring coupling constant 
g'^/N in the fermionic model is mapped into the inverse mass spin 1 massive excitation. In the same fashion, we may 
identify each term of the generalization to higher-order derivatives studied below, with terms following from the 1 /M 
expansion of the fermionic determinant. This is however beyond the scope of the present work. 



III. GENERALIZATION TO HIGHER-ORDER DERIVATIVES 

In this Section we use the duality procedure which connects the SD and MCS models, and the MCS-Proca and 
MCS-Podolsky models, in order to extend the formalism to generalized models involving higher-order derivatives. 
Evidently, the presence of higher-order derivatives introduces longer distances effects and is of interest to String 
Theory, and to investigations involving phase transition. We split the subject into two subsections, the first exploring 
the relevant classical issues, and the second dealing with duality at the quantum level. 

A. Duality at the clasical level 

Let us rescale fields and coordinates as A ^ m^^^A and x — + m^^x, in order to work with dimcnsionlcss quantities. 
We rewrite A and x as A and x again, and we define the general field A^"'' through the recursive relation 

^(n)M ^ ef'^Pd^A^^-^^ (19) 

Here we use A^^-* = A^ to represent the basic field. The subscript n which identifies the field also shows the number 
of derivatives one has to perform in the basic field. We notice that in the above relation (|19|l the field A^"'' is related 
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to the field A^"'^'^ by means of a derivative, which maintains the mechanism we have identified in the former Sec. ^ 
where the dual relation between the fields involves a derivative. 
For this field we now define the general field strength, 

F^;:)^a^4")-9.4"). (20) 

With this, the Maxwell term is proportional to A'j}'' A^^^^^ , and can be generalized to A^^'' A^"-^^^ . 
We define the antisymmetric tensor g|i"^ such that G^°j = F^^,, and 

G^-J ^ di,d'Gtj'^ = d,d'Gt~'^ - d.d'G%-'^ (21) 

which is important to relate A'jT^ to A^^K We notice that G^"^ contains 2n + 1 implicit derivatives in A^. So, the 
general field is, for n positive, 

( (-l)ta'^ g1"/^"^^ for n even 

4") = i (22) 
[i(-l)^£/^ g[,"/'-'/'^ for n odd 

To include CS-like terms into the proposed generalization we first notice that the CS term can be written as 
e'^'^^A^d^Ax = Thus, it induces the general form, A^^^A^-'^'^, where j,j are non-negative integers. We 

notice that if i = j we obtain the generalized Maxwell term. Thus, we can modify the generalized Maxwell model 
appropriately, to include extended CS contributions. Furthermore, we can prove the identity 

4)^(j> ^ 4^i)a(^'+i)^ + e^"'^^^ (4"^^^A^) (23) 

The last term is a total derivative, which can be discarded since it does not change the action when this identity is 
used in the Lagrange density of the corresponding model. 

We see that for i + j even we can write, apart from a total derivative, 

40^0-)^ ^ Ai'^^A^'^^^' (24) 

As a result, in the action we do not need to split A^'ji^ A^^^^^ into generalized Maxwell and extended CS terms. We 
introduce a single term, characterized by the sum i+ j: for i + j even we get a generalized Maxwell term, and for i + j 
odd we obtain an extended CS term. For example, and reproduce the standard CS and Maxwell 

terms, respectively. 

We use the above results to introduce the model 

n 

= ^ r,A|,°UW'^ (25) 
where n is integer, and are real and dimcnsionless parameters, with ^ 0. The equation of motion is 

n 

^r,4)=0 (26) 

This equation allows writing cJ^'ylJ^-' = 0, which shows that no longitudinal mode propagates in the theory. Besides, 
we can use this condition to rewrite Eq. (|22|) in the simpler form 

{(— l)^n"5^A|j'^^ for n even 

(27) 
(-l)VnV£^^;^a''^(o)^ torn odd 



With the aim to build an Abelian gauge model, we consider the variation 5A^ — SA'"^'' = i9^A; as before, A is a 

! definition for yl^"^ to obtain 

54) =0, z^O (28) 



local infinitesimal parameter. We use this and the definition for a']^' to obtain 
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We vary the Lagrange density in H25() to obtain 

^ i^A'J^W (29) 



where the local Noether current is defined as 



K,^2j2nA^;^ (30) 



1=0 



We introduce an auxiliary vector field a^, and we couple it linearly to the Euler vector 

£(") ^ _ gM;^^ (31) 

We chose (5a^ = 9^ A = JA^. Thus 

SC["'^ = -af.SRf' (32) 

We consider 

= + roa^a^ (33) 
The procedure ends with the elimination of the auxiliary field a^. We get to 

which is the Lagrange density of the dual model. We use the Euler vector to get 



(34) 



n ^ n n 

= ^r,4°)^«- _ (35) 

8=0 1=0 i=0 

or better 

n n 

^--VT^^^^^^^'''- (36) 



'^0 i=l j=0 



which gives the dual model. We notice that the respective action is gauge- invariant. This result shows that the model 
(I25|l engenders hidden gauge invariance, in a way similar to the SD model 0- We can rewrite this result in the form 



'2n 



i=l 

where we have used 

' 1 ' 

n Ti-i for 1 ^ i 



^0 1=1 



r = < 



(38) 



1 " 

y ri r.i-i for n < i ^ 2n 

ro ^ 



l=i — n 



Before ending this Section, let us comment on issues related to the above duality investigation. The equation of 
motion that follows from the dual model H37|l is 

2n 

J2r[A^^0 (39) 

1=1 

In Eq. (|2()(l we eliminate A^^^ to see that these two equations are identical, thus confirming the dual equivalence 
between the two theories. 
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Another issue concerns finding a master theory. To get to this, let us introduce another field B^^^ such that 
Bj!^^ = B^. We use A^^^ and B^'^ to write 

n n 

= ro4"'i?(°)'^ + 2 J2 nBj^^A^'^'^ - J2 (40) 

i=l 1=1 

We vary the corresponding action with respect to B^"^ to get 

b(0) = _ly^^^WP (41) 

We use this into H4Q|) to obtain 

^ n n 

4'^ = -rEE-^^.4"^'^'^"- (42) 

which is Eq. If we vary the master action with respect to aJ^-* to get 

n n 

^r,;i?W=^r,A« (43) 

1=1 1=1 

This result allows writing 

n n 

C = ro4")i3(")^ + J2 nBj^^B^'^'^ = nB^°^B^^^ (44) 

i=l 1=0 

which reproduces Eq. H25|l . confirming that Eq. H4U|) is a master or parent theory. 

The model that we have introduced is defined by Eq. H25|l : it is controlled by the integer n which we choose to be 
n — 1,2,3.... The simplest case is n = 1, which reproduces the SD model. The next case is n = 2, which gives the 
MCS-Proca model investigated in the former SecE] The case n = 3 gives 

£(3) ^ ro4"U(o)^ + ri40)A(i)^ + r24°)A(2)A' + r34°U(3)A' 

= roA^^A'^ + rie'^-^A^d^Ax + ^rsF^.F^'' + r^e^"^ A^d^d^ F^p (45) 
and the other cases follow standardly. For n — i the dual model is 

4'^ = --EE -^-.4^^^'^"^ (46) 

^° 1=1 j=o 



We use Eqs. and to write 



^n^ = -"^e^^^^KFuX - ^ f ^2 + F^,F^^'' + \{rz + 2^ ) e^^''^F^,^PF, 



-{rl + 2rir^)d^,F^^^d^Fx, - 2'^e^'''^dP Fp^d.d'^ F^x - i^df^d" F^xd^dpF^^ (47) 
ro ro ro 

In this example, we can choose the parameters tq 7^ 0,ri,r2,r3 to find the dual theory to every model, up to the 
order n — 3. For instance, if we use Tq = 1/2, ri = —1/2, r2 = —a/2, and = 0, and if we re-introduce dimensional 
units, we obtain from (|45|l and H47|l the former expressions and © that we have investigated in Sec. HTl 

We notice that the presence of the CS and CS-like terms imposes the restriction, that the Minkowsky space must 
have (2, 1) spacetime dimensions. Thus, if we set = for i odd, we eliminate all the CS and CS-like terms, and our 
results are then valid in arbitrary (d, 1) space-time dimensions. 
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B. Quantum duality 



We now explore the issue of quantum duality l33l |3J, |35j for the models examined previously. We start with the 
Lagrange density for the master theory, given by Eq. 14U|) . We extend the model to write 

C^m\J^J) = r-oS(,"'B(o)'^ + 2B'i°^F^{A) - A'^°^F^{A) + aroB(,o)j("'^ + aFt'{A)j(-°^^ (48) 
where we have set 



i=l 



(i) 



(49) 
(50) 



to ease calculation. We see that the master theory contains two fields, thus we have added the external currents 
and J^; we see that couples with B^, while couples with A^ indirectly, through Ffj_{A). 
We use to eliminate A^. We obtain F^{A) = F^{B) + (a/2)F^(J), which we use to write 



/:(")(J, j) = + aB(")(ro/")'' + F^(J)) + ^ j(o)i^^(J) 



(51) 



where is the model given by Eq. We also use ^ to ehminate 5^°^ We get S^"^ = -(f/ro)F/,"^(A) 



j(0) _ 



(a/2)j^^\ which we use to abtain 



4')(J, j) = + aF,(A)(j(")- - /°)-) - ^j,f 



(52) 



where Cj^' is the dual model, given by Eq. (|36|l . 

We define the functional generators z'^\j,j), Z^^^^J^j) and Z^\j,j) in the usual way. They allow obtaining 



f S'Z 



2 yin) 



M 



and 



2 



1 J2^(«) 



,j=o ^(")'5j^(x)<5>(2/) 



s^z'i!;^ 



The master theory allows writing 



f 6^Z 



2 ^(n) 



Z^^ 5j^i{x)53u{y) 



and 



1 



1 5'Z 



2 7(n) 



a\l < S("'^(a;)i3(")''(j/) > 



<F>'{A){x)F''{A){x) >D 



D 



S2 yW) 



1 J2^(„) 



W^(a;)<5J.(2/) 



j=o <5J^(x)<5J.(y) ^. 

where the operator O^''' is defined according to 

We use these results to write 



2 

-a^ < i^''(A)(x)F^(^)(y) >Z3 -i^<5'^^(5(x - y) 



^ -a^ < F''{B){x)F''{B){y) > +t—0'"'Six - y) 



< Bj^Hx)B(J>^{y) >=< F,{A){x)FM){y) >d +^^S,J{x - y) 



and 



(53) 
(54) 
(55) 

(56) 

(57) 
(58) 



< F^'{B){x)F''{B){y) > -i-O^^ b{x - y) =< F^iA){x)F'' iA){y) >d 



(59) 
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Expressions H58|l and (|59|l show that the correponding Green functions are equivalent, apart from contact terms. 
They show that the duahty that we have presented in Sec. (III. A) above is preserved at the quantum level. 

We now examine the specific case considered in Sec. El which involves duality between the MCS-Proca theory 
and the generalized MCS-Podolsky model JSjl. In this case the results H58II and (|59|l become 

9 

772 

< B^{x)B,{y) >=< F^{A){x)F,{A){y) >d +i—Sf,.Six - y) (60) 

and 

< F^{B){x)F,{B){y) > (-ye^.A^^ + ^5^.0 - ^9^,9.) b{x - y) =< F,{A){x)F,{A){y) >d (61) 

where we have re-inserted the parameters used in Sec.|n] We see that for a = one recovers the result of Rcf. [s^ . 
which deals with quantum duality for the CS and MCS models. 



IV. THE CASE n^oo 



The success of the former investigations has lead us to think on extending the case containing a finite number of 
terms to the case where an infinity sequence of terms is considered. Our interest relies on the case where the infinity 
sequence of terms adds to give elementary functions. We see that the form of in H25|l suggests this particular 
generalization, to the case where 77 ^ oo. Such generalization can be written in terms of a smooth function, depending 
on the specific values of the real parameters r,; that we have introduced to define the model. We further explore this 
possibility introducing the model 



C ^ A^lFm^^^A-" 
where is a smooth function, and O is the operator 

0^"' EE --e^"'^^x 

The model is defined in terms of the expansion of the smooth function, in the form 



n=0 



where C„ are given by 



1 d"F{x) 



i\ da;" 



x=0 



(62) 
(63) 

(64) 
(65) 



and [0°]'^ = (5'^^, [O^]^^ = -^''ux^^ and [0^]1 = [O^f^p^]),, and so forth. We require that F(0) ^ 0, which implies 
that Co 0. This means that the above model starts with the Proca-like term Afj,A^, which prevents the presence of 
gauge invariance. Thus, we use the gauge embedding procedure to obtain the dual model. It has the form 



oc oc 



n=0 



or, formally. 



C 



D 



A., 



F{0) 



We illustrate this general result with the example 

C^A, 

In this case the dual theory reads 



m 



1 

l + O 



o 



[F{0)f 



A"" 



A" 



A"" 



(66) 



(67) 



(68) 



(69) 
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V. ADDING FERMIONS 



The study of fermions is motivated by the possibhty of extending the present duaUty procedure to more reaUstic 
models, which should necessarily contain fcrmionic matter fields to describe the matter contents of any realistic model. 
We add fermions with the modification 



where 



(70) 



(71) 



describes free fermions, with M being the dimensionless fermion mass parameter. To identify how the fermionic field 



interacts with the fields A^""^ we introduce the non-minimal coupling 

n 



(72) 



where are (dimensionless) coupling constants - we notice that the condition a = Q, i ^ Q leads to minimal coupling 
with the fermionic matter. With the above generic coupling, the interactions terms in Ci are given by 



i) 



where we have defined the general current 



with 



4o) = j^ = vi7^v 



To write the above expressions we have used the identity 



(73) 

(74) 
(75) 
(76) 



We now search for the dual theory. The procedure follows as in the former case. The Euler vector is modified by 
the presence of interactions; it changes to 



K, = K,-Y,eiJ'tl 



(77) 



i=0 



We introduce an auxiliary field a^, and we impose that Sa^, = = 5K. The Lagrange density varies according to 

since (5Jji"'' = 0, as one can verify straightforwardly. Like in the former case, the procedure requires 
another iteration. The final result is 



An) _ p{n) ^_Xt^K 



or, explicitly, 



where we have set 



2n 2n 
i=l i=0 

1 * 

> n Ci-j for 1 ^ z ^ n 

1 " 

> ri Bi-i for n < i < 2n 



(78) 



(79) 



(80) 



l—i—n 
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and 



— e; ei-i for < « ^ 

^ 1=0 



(81) 



1 " 

- E 



; Bi-i 



for n < i ^ 2r 



[=i—n 



In this result, we notice the presence of Thirring-like interactions 'S&l, which are fundamental to maintain the 
contents of the fcrmionic sector unchanged [s^l • To see this we examine the dynamics of the fermionic sectors in both 
theories. The fcrmionic equation of motion for the first theory is 



i=0 



To eliminate the gauge field, we notice that 



n 1 



i=0 i=0 

We restrict the coupling constants to obey ei = ari to get 



n 



(82) 



(83) 



(84) 



i=0 

Analogously, the fermionic equation for the dual theory is given by 

2n 2n 
i=l i=l 

We find the equation of motion for the gauge field in the dual model as 

2n ^ 2n 

2 



(85) 



E-^4^ = ^E^^^^ 

i=l i=l 

The restriction = ari implies that = ar^, and so we can write 

1 " 



(86) 



(87) 



i=0 



which is Eq. H84() . This result shows that the fermionic sector does not change when one goes from H70|) to the dual 
theory ifT^ . 

We can also add fcrmions in the case n — > oo. To illustrate this possibility we consider the model 

C = A^[F{Or,A-^ - A^[G{0)r,r + Cf (88) 
where G is a smooth function similar to F. We use the gauge embedding procedure to obtain the dual model 



Cn = A„ 



A" - A, 



G{0) - J^^F{0)G{0) 



- - / 



1 



F(0) 



r + Cf (89) 



We notice the presence of the generalized Thirring-like term in the dual theory. 
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VI. NONLINEAR INTERACTIONS 



We now return to the Born-Infeld [31 generalization, which is different from the Proca and Podolsky generahzations. 
The main ingredient now is nonhnearity, and so we further explore the duality procedure introduced above mixing 
higher-order derivatives and nonlinearity. Evidently, there are several different possibilities of including nonlinearity 
in the model introduced in the former Sec. IfffI but here we consider the case 



(90) 



where g{x) is nonlinear in x — This case does not include nonlinear interactions that involve derivative of 

the basic field — A'"^^ . The equation of motion is 

A?) 



where g'{x) = dg/dx. This equation allows writing 



(91) 



(92) 



which shows that there are longitudinal mode propagating, due to the presence of the nonlinear interaction. We notice 
that in the linear case [g{x) = x] we have g' — 1, which leaves no room for propagation of longitudinal modes. 

We treat the presence of nonlinearity invoking the trick used in Ref . , and then further explored in [2^ Is^ . The 
key point here is to remove the nonlinearity under the expense of introducing another field, an auxiliary scalar field 
(j). We implement this possibility with the change 
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(^(o)^(o)p) ^ /(</>) + 140)^(0)'^ 



We follow ^1 to show that 



The model is modified to 



'dx-,g'-' (- 
X \x 



1 " 

4"^ = /w + i4U(°)-+^r.4)^( 

In this case the Euler vector is given by 



i=l 



The gauge embedding method allows writing 



and so 



where fi is given by 



2ri 



(93) 



(94) 



(95) 



(96) 



(97) 



(98) 



4=1 



/ i-1 



1=2 



n r.i^i for 2 ^ i < n 

1=1 

n 

n ri-i for n < i ^ 2r 



(99) 
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We notice that in (|98|l the nonlinear behavior involve all the terms, unless the Chern-Simons one. This fact is clearer 
when we eliminate the auxiliary field (j) from the model, which is formally given by 

<^ = /"' (100) 
We illustrate the above investigations with the model 



4") =;32y^^|j^(0)^(0)M + ^,^4°)A('> (101) 



1=1 



where the nonlinear contribution is of the Born-Infeld type. We notice that the limit /3 — > oo restores the original 
model (|25f) . In this case, the function /(0) given by Eq. H94|l has the form 



/(</>) = A) (102) 



2 V ^00 



We use Eq. ifTUUI) to obtain 




2 



^n4°^AWM (103) 



The gauge embedding procedure given above allows writing 



r(n) _ o2 
'-BID ~ 



2 



P,A''°^ A(.^)t- - ^.4°^^^'^'^ (104) 



2 

i=2 



We see that in the limit /3 ^ oo the above result leads to (|37() . the dual of the model (|25|l . as expected. We also notice 
that in the case n = 1 the model IjlOll) reproduces the Born-Infcld-Chern-Simons model investigated in Ref. 28], and 
this shows that the model HlOlfl is a generalization of the model investigated in . 

VII. COMMENTS AND CONCLUSIONS 

In the present work we have investigated duality symmetry in generalized Field Theory models, involving higher- 
order derivative of the basic field — A^^"* through the recursive relation aI""* — e^'^^d^A^"^ The investigations 



started in Sec. |IT| and in Sec. IIIII and there we have generalized the self-dual model to include several higher-order 
derivative terms, and we have obtained the dual theory. We have also proposed a master model, from which one 
gets both the model and its dual partner. In this generalization, the presence of the CS and CS-like terms imposes 
the restriction that the Minkowsky space has (2, 1) spacetime dimensions. However, we can eliminate all the CS and 
CS-like terms with the restriction = for i odd; in this case our results are valid in arbitrary (d, 1) space-time 
dimensions. In Sec. IIIII we have also examined the implications of duality at the quantum level. These results are 
obtained for for n integer, finite, and in Sec. II VI we have further extended our results, considering the limit n ^ oo, 
which leads to the case involving non-polynomial functions. 

In the sequel, in Sec. we have added fermions to the system, and there we have included fermions in models 
involving A^^\ in the case where n is finite, and also for n — > oo. In Sec. IVII we have investigated the more general 
case, where nonlinear contributions involving the basic field A^^'' is present. As we have shown, in this case we first 
circumvent nonlinearity at the expense of introducing an auxiliary field, and then we proceed as before, to get to the 
dual model. After getting to the dual model, we then eliminate the auxiliary field to obtain the dual model in terms 
of the original field aIi"-* and the accompanying derivatives. 

We notice that all the generalizations we have investigated involve the Abelian vector field A^. Thus, a natural 
and direct issue concerns the use of non-Abelian fields, and their possible generalizations, along the lines of the non- 
Abelian SD model, and the Yang-Mills-Chern-Simons model [2^ . Also, it is of interest to investigate if the presence 
of bosonic matter changes the duality scenario that we have presented in Sec. Another point concerns duality in 
models of the BAF type, and their extensions to include higher-order derivatives. These and other related issues are 
under consideration, and we hope to report on them in the near future. Another point concerns models which include 
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fermions. As one knows, if there is no fermionic self-interaction, and if one integrates on the fermions, the remaining 
action will necessarily contain higher-order derivative terms, thus giving another compelling motivation to investigate 
models involving higher-order derivatives. We hope to report on these and other related issues in another work. 
We would like to thank CAPES, CNPq, PROCAD and PRONEX for partial support. 
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